Normal state of highly polarized Fermi gases: Full many-body treatment 
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We consider a single J, atom within a Fermi sea of f atoms. We elucidate by a full many-body 
analysis the quite mysterious agreement between Monte-Carlo results and approximate calculations 
taking only into account single particle-hole excitations. It results from a nearly perfect destructive 
interference of the contributions of states with more than one particle-hole pair. This is linked to 
the remarkable efficiency of the expansion in powers of hole wavevectors, the lowest order leading 
to perfect interference. Going up to two particle-hole pairs gives an essentially perfect agreement 
with known exact results. Hence our treatment amounts to an exact solution of this problem. 

PACS numbers: 05.30.Fk, 03.75.Ss, 71.10.Ca, 74.72.-h 



The last few years have seen a very strong development 
of works, both experimental and theoretical, devoted to 
the exploration of the physical properties of ultracold 
fermionic gases JJ. In particular through the existence 
of Feshbach resonance these systems provide a physical 
realization of the BEC-BCS crossover, and around the 
resonance very simple examples of strongly interacting 
fermionic systems, of high interest for example in con- 
densed matter physics. More recently much efforts have 
focused on systems where there is a possibly strong im- 
balance between the two fermionic populations present 
in these gases [ij, l2], corresponding for example to the 
two lowest energy hyperfine states of ^Li or ''"K. 

A very interesting limiting case of such imbalanced 
mixture is the case of a single fermion (say |-spin) with 
mass m^ in the presence of a Fermi sea of another, say 
I, fermion species with mass mp This kind of mixture 
has anyway to be considered in the analysis [3|, |4| of ex- 
periments on trapped gases. Of major importance for 
the analysis of the phase diagram [3|, y, l5( are the bind- 
ing energy Eb and the effective mass m* of the |-atom. 
At unitarity recent fixed-node Monte Carlo (MC) cal- 
culations [3] give {5/3)Ei,/Ef = 0.99(1) where Ep is 
the Fermi energy of the t atoms. On the other hand 
a simple T-matrix analytical calculation, which happens 
to coincide with a variational calculation [J, \^, gives 
Ef, = 0.6066£'i? which is remarkably close to the MC re- 
sult. This closeness has been confirmed very recently by 
diagrammatic MC calculations, leading to Ei, = O.QlSEp 
[31 or £'(, = 0.615ii^_F tS|j depending on the specific de- 
tailed handling. The very close proximity of the ana- 
lytical and the MC results, considered to be quite near 
the exact result, is a major puzzle. Indeed at unitarity 
the system is strongly interacting, whereas the analytical 
treatment considers only single particle-hole excitations 
(of the non-interacting Hamiltonian He). This descrip- 
tion is only appropriate for a weakly interacting system 
and should fail to a large extent at unitarity. 

In this paper we solve this puzzle and at the same time 
provide an essentially exact solution to the problem. This 



is done by considering states with any number of particle- 
hole excitations. While the ground state of the interact- 
ing Hamiltonian H has important weights coming from 
many particle-hole excitations, we show that with respect 
to the calculation of the energy of the J,-atom these states 
give contributions which decrease extremely rapidly with 
the number of particle-hole excitations. This is basically 
due to destructive interferences between contributions of 
these many particle-hole excitations states when two par- 
ticles (or holes) are exchanged. The efficiency of these 
interferences is directly linked to the key ingredient of 
our solution, namely the fact that an expansion in pow- 
ers of the holes wavevectors q^ turns out, quite surpris- 
ingly and unexpectedly, to be an excellent approximation 
scheme. To lowest order in this approximation, i.e. when 
the dependence on the q^'s is neglected, the interference 
is perfect and there is a complete decoupling between 
the states with different number of particle-hole excita- 
tions. In particular this justifies the consideration of sin- 
gle particle-hole excitations and explains why the above 
analytical result is so good. Improving the approxima- 
tion by taking more properly into account the q^ depen- 
dence provides a small coupling between the states with 
one and two particle- hole excitations, providing a small 
correction to the preceding result. In turn an even bet- 
ter handling of this dependence gives a coupling between 
states with two and three particle-hole excitations, giving 
a small change to the preceding small correction, and so 
on. In this way we obtain a series of results with strongly 
increasing accuracy, converging extremely rapidly to the 
exact result [9|. 

Since the theoretical solution of this problem has 
reached a high and secure level of agreement, which is 
quite unfrequent in many-body problems including the 
case of cold fermionic gases, it provides a useful and 
convenient benchmark for experiments. Fortunately the 
agreement with experiments seems reasonably good |lO| . 
Moreover our solution is of high general interest since 
exact solutions of many-body problems are quite unfre- 
quent, mainly for three dimensional cases. Our finding 



that a hole-wavevector expansion works so well is also 
quite interesting, and it would be worthwhile to apply it 
to more complex many-body problems. Another useful 
conclusion regarding self-consistent treatments can also 
be taken out of our solution. An attractive and standard 
way to try to improve on a siniple T-matrix approxima- 
tion, such as the one used in [6|, is to use full propaga- 
tors rather than bare ones in the T-matrix. The practical 
handling is much more complicated, but this seems much 
more satisfactory since this amounts to include an infi- 
nite set of higher order diagrams. This has been done 
for example recently in the BEC-BCS crossover for equal 
spin populations ll|. We have made the corresponding 
self-consistent treatment for the T-matrix approximation 
of Ref.y . The result at unitarity is Eb = 0.68 Ep, much 
worst than the non self-consistent one. Our solution al- 
lows to understand this failure. Indeed the exact solu- 
tion displays almost complete cancellations between di- 
agrams, whereas the self-consistent treatment destroys 
this delicate balance since it retains only some diagrams 
at a given order. Hence self-consistency does not system- 
atically lead to an improvement. 

Here we concentrate on the unitarity case. Indeed in 
the weak limit a — > 0_, an expansion in powers of a 
should be valid and the lowest order correction should 
give already quite good results. On the other side of 
unitarity, one can easily see from the details given be- 
low that our approximation should improve for increas- 
ing 1/a > 0. Hence unitarity is a kind of worst case 
for our approach. Since it attracts in addition most of 
the interest in the strong coupling regime of large a, this 
justifies that we restrict our discussion to this case. Gen- 
eralization away from unitarity is quite obvious. 

Before presenting our solution, we display very ex- 
plicitely the efficiency of the q^ expansion on the sim- 
plest case of the lowest order solution [6(]. At unitarity, 
for equal masses tti^ = mi = m, the binding energy 
Eb = —E = \E\ is given by: 
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Ck-q + Ek - Cq with ek = k /2m. 
Setting q = in E\^^ leads to the equation 2/3p = 

1 -I- \/ p/2 arctan yj p/2 for p = \E\/Ef- The solution 
p = 0.5347 is not so far from the actual result, which 
is already pretty satisfactory for such a crude approxi- 
mation. However we may quite improve on it by treat- 
ing q to lowest significant order instead of neglecting 
it completely. This means we take for the k angu- 
lar average {[\E\ -Hek-q + ek -eq]~^)k '=^ [2efc + \E\]~'^ + 
q^k'^ /{2>w?)[2ek + \E\]~ . This gives at unitarityp — 
0.5985, which is remarkable compared to the exact [6| so- 
lution p = 0.6066 of Eq.Q and the MC result (5/3)p = 
0.99(1) [5|. If we do not restrict ourselves to unitarity 



and let the scattering length vary, the corresponding so- 
lution for p turns out to be essentially undistinguishable 
on a graph from the exact numerical solution [6|. Al- 
though the situation is somewhat less favourable for un- 
equal masses, this shows very simply and explicitely the 
validity and efficiency of this (\i expansion approach. We 
apply it now to the full many-body problem. 
The Hamiltonian of our problem is: 



H = Hc + V 
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where ek = k^/2m|, £'(P) = P^/2mj, and Ck and 
cj^ are annihilation and creation operators for |-spin 
atoms while 6p and 6p are for the j-spin atom. In 
the potential energy term V, the short-ranged interac- 
tion potential provides an upper cut-off kc, which we 
will as usual let go to infinity while the coupling con- 
stant g goes to zero, keeping the scattering length finite 
in the relation mr/{2TTa) = g~^ + X] "^ im^/k'^ , where 
TO,. — m,im,il {m^ -f mi) is the reduced mass. Momen- 
tum conservation in the scattering is insured by the Kro- 
necker symbol (5kk'PP'- Physically this potential energy 
term creates (or annihilates) a single particle-hole pair 
from the Fermi sea, or merely scatters particles or holes, 
the momentum transfer being taken by the |-spin atom. 
Hence we are looking for the ground state as a general 
superposition of states with any number of particle-hole 
pairs: 
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where |0) = rife<fc '^k 1''^'^^) is the Fermi sea of f-spins. 
For simplicity, in all the terms, we have not written ex- 
plicitely the creation operator 6p of the J,-spin since its 
momentum P is just obtained to insure that the total mo- 
mentum of the particles in any term is zero. We assume 
also implicitely ki > kp and qj < kp. The coefficients 
ttkiqj (which is a shorthand for a{ki}{qj}) are naturally 
antisymmetric with respect to the exchange of any their 
arguments k^ or q^. 

Writing H\ip) — E\tp) (where we take as zero for the 
energy E the energy of the whole Fermi sea, and omit the 
average potential energy which disappears for g —f 0) and 
identifying the coefficients of specific particle-hole states 
in both terms, we find for the equations corresponding 
to the full Fermi sea and the Fermi sea with a single 
particle-hole pair: 
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where we have introduced E^'' — \E\ + E{J2i k^ — 

X)? 'ij) + S" ^ki — X)? ^qj • ^^ tti^ second equation the 
first term conies from the creation of a particle-hole pair 
in the Fermi sea, the second from the diffusion of the 
particle, the third from the diffusion of the hole and the 
last one from the annihilation of a particle-hole pair. In 
the last term, from two particle- hole pair states, there 
are four combinations for annihilation of a particle-hole 
pair (this would be n^ in the general case). The next 
equation from two particle-hole pairs states is: 

— 1 (2) 
~ 9 ^kk'qq''^kk'qq'=-akq- ak'q'+akq'+ "k'q 
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with physical origins for the various terms analogous to 
the preceding ones. One could write formally the gener- 
alization for any order, but we will not do it since this 
will not be necessary. 

In Eq.© the second term turns out to give a diver- 
gent contribution when kc —* oo and the limit is finite 
only after multiplication by g. By comparison the third 
term, where the summation is on the Fermi sea, displays 
no divergence and gives a zero contribution after multi- 
plication by g. Hence we omit it from now on. Similarly 
in Eq.® we can omit the two analogous sums over Q. 

We see that, to solve our problem, the only thing 
we need to know from Eq.® is X^Ko'^kKqQ, which is 
the only quantity entering Eq.([71) from the subspace of 
two particle-hole pairs added to the Fermi sea. Now, if 
we neglect the wavevectors q and q' in E''-^'> and take 
we can sum Eq.® over q'. Because 
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of the antisymmetry of ctKkk'Qqq' in the exchange of Q 
and q', the last term drops out. Hence, through this 
summation, Eq.® provides us with an integral equa- 
tion for the function X^kq Q^kKqQ- We see that our ap- 
proximation makes the subspaces with higher number of 
particle-hole pairs irrelevant for our problem, and we are 
completely decoupled from these higher subspaces. We 
have obtained, together with Eq.® and Eq.Q a closed 
set of equations for ao, akq and J2k.q '^kKqQ7 which can 
clearly be solved exactly by numerical means to provide 
the energy |_B|. The precise justification for our approxi- 
mation is that fc, k' > kp and q, q' < kp make indeed in 
-^kk'qq' the k and k' terms dominant and q and q' ones 
negligible in a large part of the variables space. Natu- 
rally this is an approximation, the interference is not ex- 
act but not so far from it. Hence the contribution from 



three particle-hole pairs states is not zero, but it is quite 
small which is the essential point. 

The same argument applied at the level of Eq.([7]) al- 
lows now to understand the success of the lowest or- 
der approach. If we make the approximation i?) ~ 
\E\ + i?(k) + Ek and sum Eq.(I7]) over q, we see that the 
last term is disappearing, due to the antisymmetry of 
«kKqQ with respect to the exchange of q and Q, which 
leaves us with a decoupled set of equations for ao and 
X^kq'^kq, leading to the good approximation discussed 
at the beginning. Conversely we may decide to make 
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the approximation E^^'_ ~ E^/ at some higher order n. 
This provides by the same procedure a decoupling from 
higher order subspaces resulting in a closed set of n equa- 
tions to be solved for \E\. Hence our analysis leads us to 
a cascade of successive approximations which converges 
very rapidly in a controlled way to the exact many-body 
solution, providing at the same time a full understand- 
ing of the success of the lowest order approximations and 
a practical way to obtain the exact results within the 
precision we like. 

In practice the convergence is so fast, as shown by the 
quality of the lowest order results, that we will only need 
for any practical purpose to implement it to the next 
(i.e. second) order. However, in addition to our main 
scheme of approximations going from one order to the 
next, we may within a given order introduce an addi- 
tional graduation of approximations. First we can as de- 
scribed above make the approximation E-^ 
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and sum Eq.® over q'. However we can easily improve 
on this handling in the same spirit as for the lowest or- 
der approximation. Indeed the argument is still valid 
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if we make only the approximation i?kk'« 
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which takes better into account the q dependence. Fi- 
nally since we know that the last term in Eq.® brings 
only a small correction, we may just decide to omit it, but 
treat the rest of the equation exactly without any approx- 



imation on E, 



(2) 



which takes even better account of 



kk'qq" 

the qi dependence. This handling is at the same level as 
treating Eq.([T]) exactly. In the same way it corresponds 
to a variational calculation since we have restricted the 
Hilbert space to contain at most two particle-hole pairs. 
This is the most precise level of approximation we will 
implement. 

Introducing G(k,q, q') = ^X^K^kKqq', dividing 
Eq.® by E^^^, ,, summing over k', taking the kc -^ 
00,3 —*■ limit, which makes again some terms vanish, 
we find the integral equation: 
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The bracket in the left-hand side can be calculated ana- 
lytically. In practice this integral equation is easily solved 
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FIG. 1: Reduced | atom binding energy Et and effective mass 
m* , in the case of t and J. atoms with equal masses, for the ap- 
proximations of increasing accuracy discussed in the text. Ef 
is the Fermi energy of the | atoms. Diamond: first order ap- 
proximation [g|. Circle: second order with qq' = 00. Cross: 
second order with qq' — qO. Square: second order with no 
q approximation (variational). For the effective mass we give 
coherently the last digits of our results to display clearly the 
trend, but they should not be taken for granted. 

by iteration, the already fast convergence being increased 
by acceleration procedures. The same handling allows to 
solve the set of Eq. (|7l9p with respect to akq- 

We have first checked the convergence of our theoret- 
ical scheme in the particular case m^ — > oo where the 
binding energy is known p| to be exactly p = 1/2 at 
unitarity. This is naturally a very convenient case. How- 
ever one can easily see that the situation, with respect 
to the convergence of the q^ expansion, is slightly less fa- 
vorable than for the equal mass case. This is appearant 
in the first order result [6|] p = 0.465 which, while being 
fairly good, is not so close to 1/2. A convenient feature 
of this case is that there are no angular integrations to 
perform in solving the various equations. At second or- 
der we find 0.481 when we take qq' = 00. This improves 
into 0.487 for qq' — qO. However with no approximation 
on E^^^ (i.e. the variational result) we find p = 0.498, 
which shows quite explicitely that essentially, at second 
order, we have already fully converged toward the exact 
result. The relative sensitivity in this case to the approx- 
imations on q may be physically understood as related 
to the lack of quantum smearing from the |-atom since 
it acts as a fixed impurity. Another situation where the 
exact solution is known is the ID case [g, [ij] . Here again 
we have checked 12| that the second order coincides al- 
most with the exact results, including for the calculation 
of the effective mass 
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We have then proceeded to the same calculations for 
the equal mass situation (see FigH]), which is markedly 
better in terms of convergence than the infinite mass 
case. For the qq' = 00 approximation we find 0.6125, 
the qq' = qO result is 0.6146. Finally the variational 
result is 0.6156. This displays quite clearly the very fast 
convergence of the q^ expansion. It is very likely that 
this last result is essentially exact [13| (say, the exact 
one is bounded by 0.6158). This is supported by how 
close our above result for infinite TO| is from the exact 
one. More precisely, as seen explicitely from the m| = oo 



case, each order brings a small correction to the preced- 
ing one. Since the second order result brings typically a 
10^^ correction to the first order one, we expect the third 
order correction to be at most of order 10""'. Naturally 
the precision of our result is beyond practical use. We 
display it only to support our claim that we have a full 
solution for our problem. 

Finally we have calculated the effective mass, with the 
same approach applied to a small nonzero momentum for 
the whole system. This produces only some minor prac- 
tical complications. We find m* /mi = 1.20. Our above 
precision on the binding energy gives us typically a ±0.02 
uncertainty. In agreement with our above findings, this 
result is fairly close to the first order value 1.17. We see 
that there is some significant discrepancies with some MC 
results [3|, |5| as well as conclusions from experiments [lOl . 
Though they are not major ones, they may be relevant 
for the detailed understanding of experiments |15]. 
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